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ABSTRACT 

Recent observations of young galaxies at redshifts z ~ 3 have revealed simultaneous AGN and 
starburst activity, as well as galaxy-wide superwinds. I show that there is probably a close connection 
between these phenomena by extending an earlier treatment of the Mbh — o relation (King, 2003). As 
the black hole grows, an outflow drives a shell into the surrounding gas. This stalls after a dynamical 
time at a size determined by the hole's current mass and thereafter grows on the Salpeter timescale. 
The gas trapped inside this bubble cools and forms stars and is recycled as accretion and outflow. 
The consequent high metallicity agrees with that commonly observed in AGN accretion. Once the 
hole reaches a critical mass this region attains a size such that the gas can no longer cool efficiently. 
The resulting energy-driven flow expels the remaining gas as a superwind, fixing both the Mbh — & 
relation and the total stellar bulge mass at values in good agreement with observation. Black hole 
growth thus produces starbursts and ultimately a superwind. 

Subject headings: accretion - quasars: general - galaxies: formation, nuclei - black hole physics 
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1. INTRODUCTION 

Most astronomers now believe that the centre of ev- 
ery galaxy contains a supermassive black hole. There 
is a close observational correlation (Ferrarese & Merritt, 
2000; Gebhardt et al., 2000; Tremaine et ah, 2002) be- 
tween the mass M of this hole and the velocity dispersion 
a of the host bulge. This strongly suggests a connection 
between the formation of the black hole and of the galaxy 
itself. 

There have been many attempts to explain this corre- 
lation theoretically. Some of these appeal to the ambi- 
ent conditions in the host galaxy (Adams, Graff & Rich- 
stone, 2001), accretion of collisional dark matter (Os- 
triker, 2000), or star captures by the central accretion 
disc (Miralda-Escude & Kollmeier, 2005). A very large 
class of models (Silk and Rees 1998; Haehnelt et al., 1998 
Blandford, 1999; Fabian 1999; Wyithe & Loeb, 2003 
King, 2003; Murray et al., 2005; Sazonov et al., 2005 
Robertson et al. 2005; Di Matteo et al., 2005; Begel- 
man & Nath, 2005) use the idea of feedback by a wind 
or outflow driven in some way by the accretion required 
to grow the central black hole. In this picture, the black 
hole eventually reaches a mass Mbh such that further ac- 
cretion (and thus growth of Mbh) is prevented because 
the outflow sweeps away the ambient gas. 
This is a very natural idea, given that we know that most 
of the mass of the nuclear black holes is assembled by 
luminous accretion (Soltan 1982; Yu & Tremaine, 2002). 
Moreover, the rate at which mass tries to flow in towards 
the central black hole in a galaxy is set by conditions far 
away, for example by interactions or mergers with other 
galaxies. These rates must be at least comparable with 
the instantaneous Eddington rates in order to produce 
the known black hole masses in the available time (par- 
ticularly at higher redshifts) . Since these external condi- 
tions cannot know what the current mass of the central 
black hole is, it seems very likely that super-Eddington 
conditions prevail for most of the time that the central 
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black hole mass grows, and indeed this is what studies of 
the cosmological history of the process suggest (cf Miller 
et al., 2005). 

It is therefore reasonable to consider the effect on the host 
galaxy of such super-Eddington accretion. Simple the- 
ory (King & Pounds, 2003) motivated by X-ray observa- 
tions of outflows from bright quasars (e.g. Pounds et al., 
2003a, b; Reeves et al., 2003) suggests that the outflow 
momentum flux is comparable to that of the Eddington- 
limited radiation field, i.e. 

M out vc^, (I) 
c 

where M out is the mass outflow rate and ^Edd the Ed- 
dington luminosity, while the mechanical energy flux is 



;M out v 2 



^Edd 



(2) 



2"""""" 2M out c 2 
The wind from the central black hole sweeps up the sur- 
rounding gas into a shell. The theory of stellar wind bub- 
bles (e.g. Lamers & Casinelli 1999) shows that this shell 
is bounded by an inner shock where the wind velocity is 
thermalized, and an outer shock where the surrounding 
gas is heated and compressed by the wind. These two 
regions are separated by a contact discontinuity. The 
dynamics of the feedback process are then determined 
by whether the swept-up gas can cool (momentum- 
driven flow) or not (energy-driven flow) on a timescale 
shorter than the flow time. In fact Compton cooling is 
likely to be effective (King, 2003: hereafter K03) un- 
til a very large (galaxy-scale) mass has been swept up 
(see eq. EU below). Then eq. states that the ef- 
fect of super-Eddington accretion is to inject the Ed- 
dington momentum into the surroundings. In a simple 
picture, no further accretion on to the black hole is pos- 
sible once tha shell of swept-up matter has acquired the 
escale velocity ~ a. The resulting theory then has no 
free parameter: remarkably, it leads to an M — a relation 
Mbh oc a 4 very close to the observed one (M 8 ~ cr| 00 , 
with M s = Mbh/M q , (7200 = cr/200 km s" 1 ) in both 
slope and normalization (K03). 
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Recent observations now make it possible to check the 
consequences of this simple model further. In partic- 
ular there is good reason to believe (Alexander et al, 
2005) that the growth of the black hole mass (revealed 
by AGN activity) occurs simultaneously with the forma- 
tion of stars in the spheroid. Moreover there is now sug- 
gestive evidence of very large-scale outflows from young 
galaxies at redshifts z ~ 3 (Wilman et al., 2005; see also 
Binette et al., 2000). A fuller treatment of the problem is 
therefore appropriate. I begin by considering the initial 
momentum-driven outflow in more detail. 

2. BLACK HOLE GROWTH 

We model a protogalaxy as an isothermal sphere of dark 
matter. If the gas fraction is f g — fibaryon/ ^matter — 0.16 
(Spergel et al., 2003) its density is 



M 2 
2nGr 2 



(3) 



where a is assumed constant. The gas mass inside radius 
R is 



M{R) = 4?r / prMr 
/o 



G 



(4) 



so this is the mass of the shell of swept-up matter at ra- 
dius R. Inside a certain radius R c (see below) the shocked 
wind can always cool so we have an outflow driven by the 
momentum rate Q 



M out v 



(5) 



where M ou t is the mass outflow rate and LEdd the Ed- 
dington luminosity. 

K03 simplified the problem by considering only the driv- 
ing by this force, and asking under what conditions the 
resulting outflow velocity reaches the escape value ~ a. 
To follow the process in detail we must include the resis- 
tance by gravity. The equation of motion is 



-[M(R)R] 



GM(R)[M + M(R)] 
R? 



4irR 2 pv 2 



M out v 



L 



Edd 



(6) 



where we have used the mass conservation equation 
AttR 2 pv = M ou t for the quasar wind velocity v, and then 
to simplify the rhs. Using Q we can simplify the 
gravity term, so 

2a 2 R 



d 



M(R)R 



M 



G 



L 



Edd 



Then using 011 the lhs we get finally 



At 



{RR) 



GM 
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M 



where 



(7) 



(8) 



(9) 



It is easy to show from (JSJ that for M < M a we have 
R < 0, i.e. the shell always decelerates. 
Now let us consider the growth of the black hole mass M. 
As argued above, this is likely to be super-Eddington at 
small M and set by conditions far from the hole. Thus 



any feedback ultimately limiting the growth of M can 
only occur when the bubble radius R is similarly large. 
This requires the wind outflow to be launched with at 
least the local escape velocity near the black hole. X- 
ray observations of outflows indeed suggest this (see the 
discussion in King & Pounds, 2003). For R > GM/a 2 
we may neglect the black-hole gravity term GM/R in 
JSJ). Then integrating twice gives 



R 2 



Rn + 2R$R$t 



2a< 



1 



M 

m7 



(10) 



where R = Rq at R = Rq with Rq some large radius. One 
could in principle connect Rq , Ro to the initial conditions 
near the black hole by numerical integration of the full 
equation ©, but this is unnecessary for our purposes. 
From HlOf) we see that for M < M a the shell reaches a 
maximum radius i? m ax given by 



R 2 



= 1 



R 2 



2er 2 [l — M/M a 



after a time 



in 



RoRo 



(11) 



(12) 



2(7 2 [1-M/M (T ] 

before stalling. Thus the bubble reaches a maximum size 
i? max in a dynamical time. For consistency we require 
GM/R max <C a 2 , or that the swept-up mass 



= M(R n 



G 



i?,„ ax > 2LM 



(13) 



should be rather larger than the black hole mass. 
Once the bubble reaches the mass M max (M), the ram 
pressure of the outflow cannot sweep up any additional 
mass. Thus in particular the shocked wind gas cannot 
add to the shell, but must fall back and orbit as it cools. 
Some of this gas forms stars and begins to grow the 
spheroidal stellar component. Depending on how angu- 
lar momentum is redistributed, some of the shocked wind 
gas is instead potentially available to continue refuelling 
the black hole, very possibly at a super-Eddington rate. 
Thus if M out 3> MEdd, gas may be repeatedly recycled 
as outflow while the black hole mass M and the bubble 
radius R grow on the Salpeter timescale t s — M/M^dd, 
increasing M max (M) on the same timescale. The most 
massive stars have lifetimes less than t s , so much of 
their gas, now enriched in metallicity, is returned to the 
ISM and is potentially available for feeding the hole once 
more. 

The continuous compression and recycling occurring dur- 
ing the black hole growth must promote efficient star 
formation. Hence this picture indeed explains why AGN 
activity and starbursts are simultaneous, as in observa- 
tions of SCUBA galaxies by Alexander et al. (2005). 
Moreover it accounts for the high metallicity commonly 
observed in gas accreting in AGN. If the supply of new 
gas for accretion (e.g. produced by merger events) is 
slow, star formation will deplete the gas supply to the 
point where accretion gradually shuts off. The swept- 
up gas (by now mostly turned into stars) must fall back, 
thus forming a spheroidal stellar bulge around the black 
hole. Such a galaxy would have Mbh lying below the 
M — a relation, only growing towards it again as a new 
merger event supplies more mass. 



3 



3. THE M-a RELATION 

Given an adequate mass supply (e.g. by mergers), M 
grows on the Salpeter timescale as explained above. 
However the shell radius i? max gets very large (whatever 
the values of Ro , Rq ) as M approaches the value where 
the denominator in (II IB vanishes, i.e. 



TfG 2 



(14) 



where k is the electron scattering opacity. In particu- 
lar, at such masses i? m ax reaches the radius R c at which 
Compton losses no longer cool the shocked wind effi- 
ciently (see K03). At this point the extra pressure accel- 
erates the shell so that it escapes the galaxy entirely (see 
below), thus finally shutting off accretion. 
Thus the black hole mass cannot grow beyond that given 
by l|14|). and this is the M — a relation. The proportion- 
ality constant is twice that of the simple treatment in 
K03, but still well within observational scatter. 

4. THE BULGE MASS 

If the shocked gas no longer efficiently cools then the 
driving force on the shell is now the total gas pressure 
P, which is larger than the ram pressure pv 2 appearing 
in @. The equation of motion thus becomes 



dt 



M(R)R 



GM(R)M ( R) =4irR2p 
R 



We now have to add the energy equation, allowing for 
the rate that energy is fed into the shocked gas, minus 
the rate of work on the ambient gas and against gravity: 



d_ 



AttR 3 3 
— 2 P 



11 T 

^Edd 



d /4tt , 
di ("3~ 



4/ 9 -£ (16) 



where I have used (2) to specify energy input from the 
outflow, r) ~ 0.1 is the accretion efficiency, and I have 
used ijU to simplify the gravity term GM(R)M(R) / R 2 . 
Now using i|15|) to eliminate P from Ijlfcifl . and replacing 
the gravity terms as before, we get 



77 T 

^Edd 



leading to 



11 T 

2-^Edd 



4 

dt 



2j> 2 
G 



M(R)R 



R-d 
2 R dt 



M(R)R 



(17) 



1 

2 RR 



3RRR 



2 R 



Wf g ^R (18) 



Using M = Mbh from (|14|l in ^Edd, this has a solution 
R = v e t with 

v 3 

V c = 3^ + Wv e (19) 



The assumption v e -C o leads to a contradiction (v e ~ 
0.01c 3> c), hence this equation has the approximate 
solution 



2rj(j 2 c 



1/3 



875(T2oo km s 



(20) 



Hence the remaining gas outside R c is driven off com- 
pletely. The total gas mass turned into bulge stars is 
thus 



M b = M(R C 



i.gxio^oA^ 2 ^) fo-7 9 3/2 M 



(21) 

(from K03). As shown in K03, this stellar mass repro- 
duces the approximate relation Mbh — 10 _3 M(, (Magor- 
rian et al., 1998) if f g , v, b all differ little between galax- 
ies. The M& — Mbh relation is known to be rather 
weaker than the Mbh — cr relation, suggesting that there 
is indeed some spread here. Equation 121|) would give 
Mb ~ fi 5 if f g ,v,b were strictly the same for all galax- 
ies, formally slightly inconsistent with the Faber-Jackson 
relation. This again suggest these quantities vary some- 
what between galaxies. 

5. SUPERWINDS 

We have seen that SMBH growth is likely to drive a 
bubble into the surrounding gas and initiate vigorous star 
formation. We noted above that the lifetime of massive 
stars is less than the Salpeter time, so the supernovae 
produced by these stars must drive vigorous mass loss 
from within the bubble. Moreover equation (|20|l shows 
that once the black hole mass reaches the value given 
by the M — a relation (|14fl its Eddington luminosity can 
in principle blow away all the mass beyond the cooling 
radius. It is tempting to see this as the origin of some 
observed galaxy- wide superwinds. 

In particular it is worth noting that superwind driv- 
ing by outflows from SMBH as discussed here is more 
energy-efficient than the conventional picture of super- 
nova driving. Thus in the example discussed by Wilman 
et al. (2005) the shell has mass ~ 10 n M Q and velocity 
~ 250 km s _1 and hence kinetic energy ~ 6.3 x 10 58 erg. 
Supernovae provide ~ 10 60 erg over 10 8 yr. However, 
the superwind is expected to cool from the shock tem- 
perature ~ 10 7 K to a few 10 4 K, in photoionization 
equilibrium with the metagalactic ultraviolet radiation 
field, as it sweeps up the intergalactic medium around 
the original host galaxy. It thus becomes momentum- 
driven once again. This means that the conversion of 
supernova input energy to kinetic energy is likely to be 
much lower than the required 6.3%. By contrast, the 
Eddington momentum flux L^d/c for a ~ 10 8 Mq black 
hole operating for ~ 10 8 yr would provide a total mo- 
mentum a few times 10 51 g cm s _1 , very close to the 
estimated value. 

This paper shows that SMBH growth simultaneously 
produces starbursts and strong outflows. Both are agents 
of mass loss, but the second is much more efficient than 
the first. Thus it may be worth looking for signs of ob- 
scured AGN activity where superwinds are seen and the 
accompanying starburst does not seem sufficient to drive 
them 

6. DISCUSSION 

I have studied the feedback probably responsible for the 
Mbh — cr relation in greater detail than in K03. In par- 
ticular I allow for the effects of gravity in slowing and ul- 
timately stalling the swept-up shell of matter around the 
growing black hole. The gas trapped within this stalled 
bubble is efficiently turned into stars and can be recycled 
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for accretion and outflow, enhancing its metallicity as it 
does so. The bubble radius grows initially on a dynami- 
cal time until it stalls, and then on the Salpeter timescale 
of the central black hole. Once the bubble reaches the 
radius where it can no longer efficiently cool, it carries 
away the gas further out in a fast, energy-driven flow. 
This simple picture accounts for a number of things. The 
stellar bulge mass is the gas mass within the cooling ra- 
dius, and agrees with the Magorrian relation. The M — a 
relation follows because the black hole cannot be effi- 
ciently fuelled once it is established. The simultaneous 
growth of the stellar and black-hole masses agrees with 
the observed AGN-starburst connection. The gas ac- 
creting in AGN must have high metallicity because it 



is trapped in the stalled bubble and repeatedly recycled 
through stars. Finally, the energy-driven outflow ending 
black-hole growth suggests an explanation for galaxy- 
wide superwinds where high efficiency is indicated. 
Clearly the approach taken here is extremely simplified. 
In particular, a fuller picture requires a better description 
of how ambient gas loses enough angular momentum to 
accrete on to the central black hole. 
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